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1 .  Summary 

The  primary  accomplishment  during  the  reporting  period  was  the  development 
of  nonlinear  fracture  criteria  for  cement-based  composites.  These  fracture 
criteria  accurately  predict  two  important  aspects  of  crack  growth  in  concrete, 
fiber  reinforced  concrete  and  other  cement-based  composites:  nonlinear  process 
zone  and  quasi -stable  crack  growth  beyond  the  peak  load.  A  preexisting  crack 
in  a  concrete  structure  propagates  in  two  stages:  subcritica!  crack  growth,  and 
beyond  the  peak  load,  post-critical  crack  growth.  The  second  stage  can  only  be 
observed  if  the  structure  is  loaded  in  a  displacement-controlled  mode.  As  a 
result  of  the  nonlinear  subcritica!  crack  growth,  when  the  fracture  toughness 
of  a  concrete  specimen  (for  example,  notched-beam)  is  determined  using  conven¬ 
tional  linear  elastic  fracture  mechanics,  different  values  are  observed  depend¬ 
ing  upon  the  size  of  the  specimen. 

A  size  independent  critical  stress  intensity  factor  and  critical  crack-tip 
opening  displacement  are  proposed  as  fracture  toughness  criteria  for  the  con¬ 
crete.  To  evaluate  these  two  parameters,  nonlinear  deformations  must  be  ex¬ 
tracted  from  the  total  displacement.  How  to  evaluate  these  two  parameters  from 
the  three-point  notched  beam  test  is  detailed  in  the  attached  paper. 

Notched-beam  specimens  of  three  different  sizes  and  five  different  mate¬ 
rials  composites  were  tested  to  develop  the  proposed  theoretical  model.  Reams 
were  tested  in  a  closed-loop  testing  system  using  the  crack  mouth  opening  dis¬ 
placement  as  a  feed-back  control.  Crack  growth  was  microscopically  moni¬ 
tored.  It  was  observed  that  the  nonlinear  critical  stress  intensity  factor  (as 
defined  here)  was  independent  of  the  size  of  the  specimen  and  was  constant  dur¬ 
ing  the  post-critical  crack  growth. 


Rased  on  the  above  proposed  fracture  criteria  (critical  stress  intensity 
factor  and  critical  crack  tip  opening  displacement),  the  fracture  behavior  of 


steel  fiber  reinforced  concrete  can  be  simulated  by  incorporating  the  bridging 
effect  of  steel  fibers.  The  mechanism  of  fracture  resistance  can  he  separated 
into  three  stages  as:  subcritical  crack  growth  in  matrix  and  beginning  of  fiber 
bridging  effect;  post  critical  crack  growth  in  matrix  such  that  the  net  stress 
intensity  factor  due  to  the  applied  load  and  fiber  bridging  force  remains  con¬ 
stant;  and  a  final  stage  where  the  resistance  to  crack  separation  is  provided 
exclusively  by  fibers.  A  theoretical  model  is  proposed  to  predict  the  fracture 
behavior  of  fiber  reinforced  concrete  based  on  the  aforementioned  fracture 
mechanisms.  The  proposed  model  is  detailed  in  the  attached  paper. 

In  addition  to  the  fracture  mechanics  based  model,  a  micromechanics  based 
model  was  developed  to  predict  the  tensile  stress-strain  response  of  fiber  re¬ 
inforced  composites.  This  model  incorporates  the  inelastic  behavior  of  the 
fiber-matrix  interface  and  the  softening  characteristics  of  the  matrix.  The 
model  is  sensitive  to  parameters  such  as  volume  fraction  and  aspect  ratio  and 
modulus  of  elasticity  of  fibers. 

Tor  fiber  reinforced  cement  based  composites,  the  bond  between  fiber  and 
matrix  is  relatively  weak  and  hence  it  critically  influences  the  performance  of 
the  comDosite.  The  bond  properties  are  generally  determined  from  a  single 
fiber  pullout  test.  To  determine  bond  strength  from  the  observed  maximum  pull¬ 
out  load,  it  is  usually  assumed  that  the  interfacial  shear  stresses  are  uni¬ 
formly  distributed  along  the  embedded  length  of  the  fiber.  An  alternate  inter¬ 
pretation  of  the  pullout  test  based  on  fracture  mechanics  concept  was  developed 
in  this  study.  In  this  model,  the  interface  is  characteri zed  by  its  specific 
work  of  pullout  fracture  rather  than  shear  strength. 
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2.  Research  Objective 


To  develop  rational  methods  to  predict  fracture  toughness  of  fiber  rein¬ 
forced  concrete  is  the  primary  goal  of  this  investigation.  It  is  known  that 
the  addition  of  randomly  distributed,  short,  steel  fibers  substantially  en¬ 
hances  the  crack  propagation  resistance  of  concrete.  How  to  experimentally 
evaluate  this  critical  property  of  fiber  reinforced  concrete  and  how  to  ration¬ 
ally  predict  it  are  the  objectives  of  this  research. 
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ABSTRACT 

—A 

A  fracture  mechanics  based  theoretical  model  is  presented  to  predict  the 
crack  propagation  resistance  of  fiber  reinforced  cement  based  composites. 
Mode  I  crack  propagation  and  steel  fibers  are  treated  in  the  proposed  model . 
The  mechanism  of  fracture  resistance  for  FRC  can  be  separated  as:  subcritical 
crack  growth  in  matrix  and  beginning  of  fiber  bridging  effect;  post  critical 
crack  growth  in  matrix  such  that  the  net  stress  intensity  factor  due  to  the 
applied  load  and  the  fiber  bridging  closing  stresses  remain  constant;  and  a 
final  stage  where  the  resistance  to  crack  separation  is  provided  exclusively 
by  fibers.  The  response  of  FRC  during  all  these  stages  was  successfully 
predicted  from  the  knowledge  of  matrix  fracture  properties  and  the  pull-out 
load  vs.  slip  relationship  of  single  fiber.  The  model  was  verified  with  the 
results  of  experiments  conducted  on  notched  beams  reported  here  as  well  as  by 
other  researchers.  Beams  were  loaded  in  a  closed-loop  testing  machine  so  as 
to  maintain  a  constant  rate  of  crack  mouth  opening  displacement. 

Introduction 


Research  conducted  during  the  last  twenty  years  has  shown  that  the  addi¬ 
tion  of  fibers  significantly  improves  penetration,  scabbing  and  fragmentation 
resistance  of  concrete.  The  possible  applications  of  fiber  reinforced  con¬ 
crete  (FRC)  include  explosion  and  shock  resistant  protective  structures.  Even 
though  the  enhanced  "cracking  resistance"  is  the  most  important  attribute  of 
FRC,  there  are  no  rational  methods  of  measuring  or  predicting  this  important 


material  property. 

For  fiber  reinforced  cement  based  composites,  the  principal  beneficial 


effects  of  fibers  accrue  after  the  cracking  of  matrix  has  occurred.  For  loads 
beyond  which  the  matrix  has  initially  cracked,  the  further  crack  extension  and 
opening  is  resisted  by  bridging  of  fibers  across  the  crack.  Many  investiga¬ 
tors  have  used  fracture  mechanics  concepts,  [1,3,5,6,11,13,18,19]  to  incorpor¬ 
ate  the  effects  of  fibers  bridging. 

In  this  paper  a  fracture  mechanics  based  theoretical  model  is  presented. 
To  aid  in  development  of  the  proposed  theoretical  model,  experiments  were  con¬ 
ducted  on  unreinforced  and  steel  fiber  reinforced  notched-beam  specimens  of 
various  sizes.  Beams  were  loaded  in  a  closed-loop  testing  machine  so  as  to 
maintain  a  constant  rate  of  crack  mouth  opening  displacement  (CMOD). 

PROPOSED  MODEL  FOR  CRACK  PROPAGATION 

The  proposed  model  can  be  described  by  considering  the  load  vs  crack 

mouth  opening  displacement  relationship  of  a  notched  beam  shown  in  Fig.  1. 

The  critical  load  carrying  capacity  of  the  unreinforced  matrix  is  assumed  to 

S  6 

reach  when  the  stress  intensity  factor  K.  is  attained.  This  K,  is  not  the 

Ic  Ic 

conventional  stress  intensity  factor  because  of  the  nonlinear  crack  growth 
that  preceeds  the  critical  load  (Fig.  1).  To  calculate  the  critical  stress 
intensity  factor,  the  inelastic  displacement  are  separated  from  the  total  re¬ 
sponse  and  it  is  assumed  that  the  effective  slow  crack  growth  (also  often 
termed  process  zone)  (*ec)  is  related  to  the  critical  value  of  the  crack  tip 

C 

opening  displacement  (CT0Dc).  The  method  of  calculating  Kj  and  CTf>nc  is 
described  later. 

Beyond  the  critical  load  if  the  unreinforced  structure  is  loaded  at  a 

moderate  rate  in  the  displacement  control  model,  then  the  crack  propagates  in 

a  "steady  state"  manner  such  that  the  stress  intensity  factor  continues  to  be 
S 

equal  to  Kj  as  depicted  in  Fig.  1. 

For  fiber  reinforced  composites,  the  crack  propagation  in  the  matrix  is 
assumed  to  be  governed  by  the  same  criteria  as  those  for  unreinforced  ma- 


n 
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trix.  However,  the  fiber-bridging  effect  must  be  included  in  calculating  the 
imposed  stress  intensity  factor  (Kj).  The  total  load  P  acting  on  the  com¬ 
posite  can  be  divided  into  three  parts  (Fig.  2): 

P  =  PM  +  PjJ  +  P[  (1) 

M  f 

where  P  is  the  contribution  due  to  matrix  and  is  related  to  Kj,  P^  is  re¬ 
lated  to  and  accounts  for  the  singularity  effect  due  to  fiber  bridging, 
and  P^  satisfies  global  equilibrium  due  to  fiber  bridging  forces  (Fig.  2). 

Note  that  the  value  of  P  for  FRC  does  not  necessarily  attain  a  maximum 
when  Kj  just  reaches  Kjc.  Depending  upon  the  volume  of  fibers,  the  maximum 
load  for  FRC  may  occur  for  a  larger  crack  length  than  that  corresponding  to 
the  peak  load  of  the  unreinforced  matrix  (Fig.  1).  Beyond  this  load,  the 
cracks  in  FRC  continue  to  extend  while  maintaining  =  Kjc  ,  the  correspond¬ 
ing  value  of  P  can  be  determined  using  the  proposed  model  (Eq.  1). 

When  the  crack  mouth  opening  displacement  (CMOD  in  Fig.  1)  becomes  very 
large,  the  resistance  offered  by  matrix  becomes  negligible  and  eventually  the 
stress  intensity  factor  (Kj )  becomes  zero.  Further  crack  separation  is  now 
,.^,nly  resisted  by  fibers.  At  this  stage  the  load  (P)  and  the  corresponding 
CMOD  can  be  calculated  from  only  the  global  equilibrium  consideration.  That 


P  =  P1 


The  beginning  of  this  stage  is  assumed  to  be  reached  when  the  CMOD  of  FRC 
reaches  CMOD  for  unreinforced  matrix  at  the  point  where  the  load  in  the  post¬ 
peak  region  is  5%  of  the  maximum  load  (Fig.  1).  At  thi s  CMOD  it  was  assumed 
that  Kj  suddenly  becomes  zero  instead  of  its  prior  value  of  K*  .  In  reality 
the  stress  concentration  value  gradually  diminishes:  this  transition  zone  was 
not  considered  in  the  proposed  model . 


EXPERIMENTAL  DETERMINATION  OF  CRITICAL  STRESS  INTENSITY  FACTOR  (K*  ) 

AND  CRITICAL  CRACK  TIP  OPENING  DISPLACEMENT  (CTODc) 

The  proposed  model  is  applicable  in  general  to  specimens  or  structures  of 
any  geometry.  To  experimentally  verify  and  numerically  demonstrate  the  model, 
notched  beam  specimens  which  were  simply  supported  and  centrally  loaded  (the 
so-called  three-point  bend  specimen)  were  tested.  To  demonstrate  the  validity 
of  Kjc  and  CT0Dc  as  valid  material  parameters,  unreinforced  specimens  of 
varying  sizes  and  made  with  cement  paste,  mortar  and  concrete  were  tested 
[7,8],  Here,  only  the  results  of  one  series  of  mortar  specimens  which  had 
somewhat  similar  composition  as  the  fiber  reinforced  mortar  specimens  to  be 
described  later  are  reported. 

Dimensions  of  specimens  are  reported  in  Table  1.  For  each  specimen  size, 
some  specimens  were  loaded  monotonically  while  some  were  loaded  and  unloaded 
repeatedly.  The  load  was  applied  so  as  to  maintain  a  constant  rate  of 
increase  of  CMOO. 

As  a  result  of  the  nonlinear  slow  crack  growth  prior  to  the  peak  load  one 
cannot  use  the  initial  notch  length  (aQ)  to  calculate  the  critical  stress  in¬ 
tensity  factor  from  the  measured  maximum  load  and  using  LEFM.  An  effective 

crack  length  (a  =  a.  +  l  )  was  defined  such  that  the  measured  elastic  crack 

0  ec 

0 

mouth  opening  displacement  (CMQD  )  and  that  calculated  using  this  length, 
► 

measured  load  and  using  LEFM  were  equal.  Once  the  effective  crack  length  is 

determined  (an  iterative  procedure  was  employed)  then  the  values  of  and 

critical  crack  tip  opening  displacement  [Fig.  3a]  (CTOD  )  can  be  calculated 

using  LEFM.  The  linear  elastic  fracture  mechanics  (LEFM)  based  formul  i  used 

to  relate  the  measured  load  (P_av),  CM0De,  CTOD  ,  and  kf  for  the  notched 

max  c  i  c 

beam  geometry  used  here  are  given  in  the  Appendix  [9,  16]. 

The  values  of  Kj  and  CTOD^  for  three  different  geometries  of  mortar  are 

shown  in  Table  1  along  with  some  other  relevant  information.  It  can  be  ob- 
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served  that  the  value  of  Kjc  calculated  as  proposed  is  essentially  inde¬ 
pendent  of  the  geometry  of  the  beam  specimens.  The  size  independency  of  the 
proposed  Kjc  was  also  demonstrated  for  beam  specimens  of  concrete  and  cement 
paste,  tested  by  the  authors  as  well  as  by  other  investigators  [7-9]. 

The  values  of  the  critical  crack  tip  opening  displacement  for  mortar 
specimens  (CT0Dc)  are  shown  in  Table  1.  The  values  of  the  critical  crack  tip 
opening  displacement  were  not  consistently  the  same  for  different  sizes  of 
beams  for  all  mix-proportions.  However,  it  was  assumed  that  they  can  be  con¬ 
sidered  independent  of  the  dimension  of  the  beams. 

It  should  be  noted  that  the  critical  strain  energy  release  rate  (Gjc) 
s  2 

is  not  equal  to  (Kjc)  /E 1  because  of  the  nonlinear  energy  absorbed  at  crack 

front  (in  the  process  zone)  as  well  as  energy  absorbed  in  the  noncritical 

region.  This  can  be  seen  when  comparing  g5  value  in  Table  1  with  the 

lc 

measured  value  of  Gp  which  was  the  total  measured  area  under  the  entire  load- 
deflection  curve.  This  additional  nonlinear  energy  was  shown  to  depend  on  the 
geometry  of  the  specimen  [7].  For  fiber  reinforced  composites,  since  the  con¬ 
tribution  of  the  matrix  fracture  energy  is  small  compared  to  that  provided  by 
fibers  pulling  out,  it  was  assumed  that  the  total  fracture  energy  of  matrix 
can  be  obtained  by  simply  multiplying  Gjc  by  a  constant  C.  The  value  of  this 
constant  obtained  from  the  ratio  Gp/Gjc  is  shown  in  Table  1. 

RESISTANCE  PRODUCED  BY  FIBERS 

The  resistance  offered  by  fibers  (only  steel  fibers  are  considered  here) 
depends  primarily  on  the  interfacial  bond  between  the  fibers  and  the  matrix. 
This  is  because  of  the  short  length  (of  the  order  of  25mm)  of  the  fibers  and 
rather  weak  zone  in  the  matrix  that  is  observed  surrounding  the  fibers  [4]. 
It  is  assumed  that  the  fiber-bridging  forces  can  be  calculated  from  the  pull¬ 
out  test  results  of  single,  aligned  fiber.  Load-slip  relationship  of  the  two 
types  of  steel  fibers  observed  from  pull-out  tests  are  shown  in  Fig.  4.  The 
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result  shown  for  Fig.  4 a  refers  to  the  fiber  with  the  hooks  at  the  ends  while 
Fig.  4b  is  for  straight  fiber.  It  can  be  seen  that  for  both  types  of  fibers 
the  initial  debonding  region  is  followed  by  fibers  pulling  out  of  the  ma¬ 
trix.  When  the  total  observed  slip  was  equal  to  the  embedment  length  (which 
was  half  the  length  of  the  fiber)  the  load  drops  to  zero.  For  both  types  of 
fibers,  the  initial  debonding  can  be  approximated  by  a  vertical  straight  line 
up  to  the  maximum  load.  With  this  approximation,  the  load-slip  curve  can  be 
represented  by: 


£M  =  n 
a  L 

max 


w 

w 

max 


0  <  w  <  w 

max 


(3) 


where  a  is  the  maximum  load  divided  by  the  total  cracked  area,  w  =  slip  or 

luaX 

crack  opening  displacement,  w  =  slip  at  zero  load  which  can  be  assumed  to 

max 

be  half  the  length  of  the  fiber,  and  o(w)  is  the  load  at  any  given  slip  di¬ 
vided  by  the  cracked  area  and  is  termed  fiber  pull-out  stress,  or  the  fiber 
bridging  closing  pressure  and  m  is  constant  which  depends  on  the  type  of 
fibers  and  was  assumed  to  be  2  for  the  straight  fibers  used  in  this  investiga¬ 
tion  [2,10,12,15]. 

The  maximum  pull-out  strength,  will  depend  on  the  number  of  fibers 

crossing  the  cracked  surfaces,  that  is,  on  the  volume  fraction  of  fibers.  In 
actual  composites,  fibers  are  not  aligned  but  have  random  spatial  distribu¬ 
tion;  this  will  influence  the  value  of  a  [12].  To  account  for  the 

m3  x 

spatial  distribution  of  fibers  an  effective  volume  fraction  of  fibers  (vef) 

was  used  in  this  study.  Based  on  single  fiber  pull-out  tests  of  straight 

steel  fibers,  the  following  value  for  a  was  used: 

max 


o  =  240  v  -  psi 
max  ef 

=  1.655  v  ,  MPa 
et 


(4) 


ef 


The  determination  of 


was  empirical  and  is  described  later. 


CALCULATION  OF  LOAD-CMOO  CURVES  FOR  FIBER  REINFORCED  COMPOSITES 

From  the  knowledge  of  Kjc  ,  CTOD^  and  the  fiber  closing  pressure  vs.  slip 
relationship  determined  from  the  pull-out  test  one  can  calculate  the  entire 
load  vs.  CMOD  and  load  vs.  deflection  relationships.  The  detailed  equations 
and  procedures  are  given  in  the  Appendix. 

The  procedure  involves  calculating  the  applied  load  and  the  corresponding 
CMOD  for  a  given  value  of  the  effective  crack  extension  (l  )  .  The  procedures 
outlined  below  are  repeated  for  different  increments  of  i  to  obtain  the 
entire  P  -  CMOD  response. 

1.  Calculation  of  critical  effective  crack  extension  (*,  ) 

v  ec 

From  the  given  material  parameters,  Kjc  and  CT0Dc>  one  can  calcu¬ 
late  the  critical  effective  crack  extension  (i  )  for  the  unreinforced  matrix 

ec 

using  LEFM  based  formulae  by  solving  Eqs.  A1 ,  A2  and  A3  for  the  notched-beam 
geometry  used.  It  is  assumed  that  the  value  of  remains  unchanged  when 

fibers  are  used. 

2.  Calculation  of  net  stress  intensity  factor  (K ^ )  for  a  given  value 
of  i 

e 

If  a  given  value  of  *  is  less  than  i  then  KT  is  less 

e  ec  1 

than  K^c  ;  and  if  >  egc  then  Kj  =  Kjc  .  For  the  precritical  crack  growth 
the  relationship  between  and  Kj  was  assumed  as  described  in  Eq.  A4  (Ap¬ 
pendix). 

3.  Calculation  of  CM0De  and  PM 

For  the  given  vlaue  of  4  and  Kj ,  one  can  calculate  PM  and  CM0De 
using  Eqs.  A1  and  A2. 

4.  Calculations  of  CMOD  from  CM0De 

Note  that  the  total  crack  mouth  opening  displacement  (CMOD) 
includes  both  the  elastic  and  the  inelastic  components.  An  empirical  relation 
was  derived  from  the  load-unload  tests  on  unreinforced  matrix.  This 


relationship  is  given  in  Appendix  (Eqs.  A5  to  A7). 

5.  Calculation  of  fiber-bridging  closing  pressure 

For  a  given  value  of  CM0I3,  the  closing  pressure  can  be  calculated 
using  Eqs.  3  and  4  and  from  the  knowledge  of  the  crack  profile.  An  assumption 
of  a  straight  crack  profile  was  found  to  be  sufficiently  accurate. 

6.  Calculation  of  Kj 

To  calculate  the  stress  intensity  factor  due  to  the  fiber  bridging 
closing  pressure  (K^)  the  procedures  described  in  Appendix  (d)  were  used. 

7.  Calculation  of  p£ 

The  value  of  the  load  can  be  calculated  from  the  known  value 
of  Kj  and  the  LEFM  based  equation  (Eq.  A2). 

8.  Calculation  of  Pf 

s 

The  value  of  P$  is  determined  using  global  equilibrium  condition 
as  illustrated  in  Fig.  A3  and  detailed  in  Appendix  (e). 

9.  Calculation  of  total  applied  load 

The  total  applied  load  corresponding  to  the  given  value  of  is 

u  f  f 

the  sum  of  P  ,  Pk  and  Pg  . 

CALCULATIONS  OF  LOAD-DEFLECTION  RESPONSE  OF  FRC 
To  relate  the  loads  with  the  load-point  deflection,  the  concept  of  global 
energy  balance  was  used.  The  total  strain  energy  release  rate  (termed  GR)  for 
the  critical  section  can  be  derived  as: 


2 

K  CTOO 

S  =  dl  (W  =  We  +  V  *  C  E^~  +  /  a(w)  dw  (5) 

where  and  Wp  are  the  elastic  and  inelastic  energies  consumed  during  the 
formation  of  new  crack,  Kj  =  net  stress  intensity  factor  due  to  applied  load 
and  closing  pressure,  CTOD  =  crack  opening  displacement  at  the  original  crack 


tip,  (Fig.  3),  and  C  =  energy  correction  factor  for  the  plain  matrix  as  indi¬ 
cated  earlier.  From  Eq.  5  one  can  calaulate  GR  for  a  given  i  (and  thus 
given  Kj,  CMOD  and  CTOD  (Eq.  A3  ))  and  obtain  the  so-called  crack-resi stant 
curve  or  R-curve  (Fig.  5). 

For  a  certain  amount  of  increment  of  effective  crack  extension,  the  total 
energy  absorbed  determined  from  GR  vs.  t  is  equated  to  that  obtained  from 
the  load  vs.  load-point  deflection  curve  (shaded  areas  in  Fig.  5).  From  this 
equality  and  knowing  the  R-curve,  load-point  deflections  can  be  determined. 
It  was  assumed  that  unloading  was  elastic  prior  to  peak,  whereas  for  load- 
deflection  curve  beyond  the  peak,  it  was  assumd  that  the  elastic  deflection 
after  unloading  remained  constant  and  was  equal  to  that  at  the  peak  (Fig.  5). 

TEST  PROGRAM 

Three  point  bend  test  was  used  to  verify  the  validity  of  the  proposed 
model.  Fiber  reinforced  beams  with  dimensions  11  in.  (280mm)  x  3  in.  (76mm)  x 
0.75  in.  (19.1mm)  (span  x  depth  x  thickness)  and  different  fiber  volume  frac¬ 
tions  (ranging  from  0%  to  1.5%)  were  prepared.  One  inch  long  (25.4mm)  brass 
coated  smooth  steel  fibers  with  0.016  in.  (0.4mm)  diameter  were  used.  Four 
different  series  with  the  same  mortar  matrix  were  cast.  The  mix  proportions 
are  listed  in  Table  2.  All  the  notches  were  saw-cut  and  were  one  inch 
(25.4mm)  deep.  These  beams  were  tested  after  curing  for  approximately  90  days 
in  an  environment  with  96%  relative  humidity  and  80°F  (26.7°C).  All  the  beams 
were  tested  in  a  servo-control! ed  closed-loop  MTS  testing  machine.  CMOD  was 
used  as  the  feedback  signal  of  the  control  system.  The  beams  were  loaded  such 

that  the  rate  of  increase  of  CMOD  was  maintained  at  2.5  x  10"4  in  (6.4  x  10" 
mm)  per  second.  There  were  6  beams  for  each  series.  A  total  of  24  beams 
were  tested.  At  least  two  beams  were  loaded  cyclically  and  the  rest  were 

loaded  monotonical ly. 
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DISCUSSION  OF  TEST  RESULTS 

The  material  properties  such  as  Kjc>  CTOD ,  a,  3,  and  C  of  the  unrein¬ 
forced  matrix  were  directly  calculated  from  the  experimental  results  of  FO 
series  and  are  775.4  psi  /T n  (0.852  MNm”^) ,  0.00025  in.  (6.35  x  10~3mm),  7, 
1.3,  and  1.25  respectively.  The  Young's  modulus  for  this  series  was  5.87  x 
106  psi  and  was  found  to  be  about  the  same  as  other  series. 

To  account  for  distribution  of  fibers  from  section  to  section  of  the 
beam,  from  beam  to  beam  and  to  include  spatial  distribution  at  a  given  section 
an  effective  volume  fraction  (vg^)  rather  than  global  volume  fraction  (vf)  was 
used  in  the  theoretical  analysis.  It  was  mentioned  earlier  that  for  very 
large  CM0D,  the  matrix  contribution  can  be  disregarded,  Kj  *  0,  and  P  =  p£. 
For  each  beam,  a  value  of  P^  was  selected  for  a  very  large  value  of  CM0D  (see 
Fig.  1).  For  this  point,  since  the  crack  had  visibly  traversed  almost  the  en¬ 
tire  depth  of  the  beam,  it  was  assumed  that  iQ  -  .99  (b  -  a^).  For  this  iQ, 

Vef  was  determined  (by  iteration  process)  using  Eqs.  3,  A9  and  Ain  (see  Fig. 
A3). 

The  ranges  of  effective  fiber  volume  fractions  (vef)  <*re  listed  in  Table 
2  in  comparison  with  the  actual  fiber  volume  fraction  (vf).  It  can  be  seen 

that  vef  ranged  about  0.5  v^  to  1.5  v^  which  seems  reasonable.  The  experimen¬ 

tal  results  of  load-CMOD  curves  for  beams  made  with  different  fiber  volume 
fractions  (including  unreinforced  matrix)  are  plotted  in  Fig.  6  and  compared 
with  the  theoretical  prediction.  The  theoetical  prediction  is  judged  to  be 
quite  satisfactory.  Good  agreement  was  also  found  between  the  theoretical 
prediction  and  experimental  results  of  load-deflection  curves  (Fig.  7).  Fig. 
8  shows  a  plot  of  the  theoretically  calculated  relative  peak  load  values  vs. 
effective  fiber  volume  fractions.  The  strength  of  FRC  beams  with  effective 
fiber  volume  fraction  of  2.5%  is  about  twice  the  strength  of  unreinforced 
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matrix.  The  calculated  values  of  energy  release  rate  [GR  >val ues  at  l 
equals  1.8  in,  (45.7mm)  (a/b  =  0.933)]  are  also  plotted  in  Fig.  8  for  differ¬ 


ent  values  of  vgf.  It  can  be  seen  that  the  energy  absorption  ability  for 
beams  with  v  ^  =  2.5%  is  about  30  times  that  of  unreinforced  matrix.  In  com¬ 
parison  to  the  improvement  of  energy  absorption,  the  strength  improvement  due 
to  addition  of  fibers  is  less  significant.  This  was  also  shown  by  Shah  and 
Rangan  [14], 

COMPARISON  WITH  OTHER  EXPERIMENTS 

Load-CMOD  curves  reported  by  Velazco,  et  al.  [17]  were  also  analyzed. 
The  material  properties  used  to  predict  the  composite  behavior  were  calculated 
from  the  reported  data  and  are  given  in  Fig.  9.  Experimental  results  reported 
by  Velazco,  et  al .  [17]  are  also  plotted  in  Fig.  9.  The  theoretical  predic¬ 
tions  are  in  good  agreement  with  the  experimental  results  for  all  different 
effective  fiber  volume  fractions. 


CONCLUSIONS 

1.  A  Fracture  Model  is  proposed  for  fiber  reinforced  concrete.  This  model  is 

based  on  the  matrix  fracture  properties  that  the  crack  propagation  in  the 
matrix  can  be  described  by  two  parameters:  critical  stress  intensity 

factor  and  critical  crack  tip  opening  displacement.  Because  of  the 
nonlinear  slow  crack  growth,  two  fracture  parameters  are  required.  The 
stress  intensity  factor  is  calculated  at  the  tip  of  the  effective  crack 
rather  than  the  initial  crack  length. 

2.  Effect  of  fibers  is  to  reduce  the  stress  intensity  factor  at  the  tip  of 
the  effective  crack  and  to  provide  additional  energy  due  to  debonding. 


17 


i  ww » •v*'v  "V  «■«  •' .  •*. 


These  effects  can  be  incorporated  in  the  proposed  model  if  the  pull-out 
load-slip  relationship  of  a  single  fiber  is  known. 


3.  The  model  predicted  load  vs.  deflection  and  load  vs.  crack  mouth  opening 
displacement  relationship  compared  favorably  with  the  experiments  des¬ 
cribed  here  as  well  as  those  reported  by  others. 

4.  In  this  investigation  only  steel  fibers  and  mode  I  crack  propagation  are 
considered. 
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Appendix 


LEFM  based  formulae  of  3-PB  tests 

For  the  notched  beam  geometry  used  in  present  study  [Fig.  Al],  the  CM0De 
can  be  calculated  as  follows  [9,  16], 

CMOO6  =-p“V1£),  a  -  a  +  *  (Al ) 

b  cE 1  10  o  e 


vt(k)  =  0.76  -  2.28A  +  3. 87A 2  -  2.04A3  +  —  6~,  A  -  ~ 
1  b  (1-A)2  b 


in  which  E'  =  - x  ,  E  =  Young's  modulus  and  \>  =  Poisson's  ratio.  The 

lV 

stress  intensity  factor  (Kj)  can  be  expressed  as: 


*i  '  tP-  Fi  (f>  ^ 

b  t 


(A2) 


The  mathematical  expression  of  crack  tip  opening  displacement  (CTO D)  can 
be  stated  as  [9]: 

a.  a. 


CTOD  =  CM0De[ (1  -  -V  +  (-1.149  £  +  1.081)(-^-  (~)Z)]  (A3) 


0.2. 


The  above  equations  are  valid  for  -g  =  4  and  gives  less  than  1%  error. 
Relationship  of  and  Kj  at . precritical  crack  growth  stage 
At  precritical  crack  growth  stage  the  effective  crack  extension  (iQ)  and 
stress  intensity  factor  (Kj)  of  unreinforced  matrix  can  be  empirically 
related  as: 

‘e  =  °'5  KIc  <KI  ‘Klc 

*Ic 


(A4) 


*e  -  0,  Kj  <  0.5  K 


where  £gc  =  critical  effective  crack  extension,  n  =  a  constant  which 


( 

i 

depends  on  the  ratio  of  a0/b  and  can  be  determined  by  curve  fitting  the  j 

I 

results  of  P-CMOD  curve  of  unreinforced  mortar.  It  was  found  that  n  =  3  1 

I 

for  the  present  study.  j 

c.  Relationship  between  total  crack  mouth  opening  displacement  (CMOD)  and  ' 

the  elastic  crack  mouth  opening  displacement  (CM0De)  ( 

I 

It  was  observer  from  experimental  results  [8]  that  the  total  crack  mouth 
opening  displacement  (CMOD)  and  elastic  crack  mouth  opening  displacement 

l 

(CM0De)  can  be  expressed  as:  ; 


CMOD  =  CM0De  ,  0.5  >  Kj 

(A5) 

CMOD  =-^-CM0De,  0.5  <  Kj  <K^ 

(A6) 

CMOD  =  CM0De(-^j-)  +  CMOD  (•—-),  KT  =  KTS 

XB-1  '  ma  x v  aB -a '  I  Ic 

(A7) 

where  a  and  B  are  material  constants  and  CM0Dmax  is  CMOD  value  at  peak 
load  for  unreinforced  matrix.  It  should  be  noted  that  a  is  indeed 
function  of  Kj  and  CTOD.  For  simplicity  a  is  assumed  to  be  constant  in 
present  study.  Values  of  a  and  B  are  reported  in  Table  1. 


d.  Calculation  of  K 


The  stress  intensity  factor  of  a  crack  in  an  infinite  strip  of  unit 
thickness  subjected  to  a  unit  point  load  (Fig.  A2)  was  used  as  Green's 
function  to  calculate  the  stress  intensity  factor  due  to  fiber  bridging 
force  (K^)  .  This  Kj  value  can  be  expressed  as  [16]: 


Kt  =  /,  2o(X)F(X,Y)//¥a  aX 

l  aQ/a 


(A8) 


where 


X  = 


Y  = 


22 


/I  -  x' 


.36  -  5.28 


+  0.83  -  1.76  X}  {1  -  (1  -X )Y } 


a(X)  is  calculated  using  the  known  CMOD  value,  crack  profile  and  Eq.  3. 
In  the  above  equation,  _a  is  given;  therefore,  X  is  the  only  variable. 
Because  of  the  singularity  problem  at  the  integral  boundary,  the 
integration  was  achieved  using  Gauss-Chebyshev  quadrature. 

Calculation  of  P^ 

The  value  of  P^  is  determined  using  global  equilibrium  as  illustrated  in 
Fig.  A3.  The  total  tension  force  at  the  cracked  cross-section  of  the 
beam  is  given  by: 

l  l 

FT  =  ;0e  to(x)  dx  +  r*  tac(x)  dx  (aq) 

e 

where  o(X)  is  the  fiber  bridging  closing  pressure  expressed  in  terms  of 
beam  depth,  is  the  stress  in  the  uncracked  portion  of  the  cross- 
section  which  is  assumed  to  be  linearly  related  to  the  beam  depth, 
and  locates  the  neutral  axis  and  can  be  determined  by  equating  the 
tensile  and  the  compressive  forces  on  the  cross-section.  The  value 
of  P^  is  determined  from: 


PI  -  “PVV  !  <b  -  Vao»]  7 


(A10) 


where  £  determines  the  location  of  the  centroid  of  FT  and  can  be 
P  1 

evaluated  from  the  shape  of  the  closing  pressure  distribution. 


TABLE  2  -  MIX-PROPORTION  OF  FIBER  SERIES 


(a)  Effective  Griffith  crack 


(b)  Crack  with  closing  pressure 


Energy  Consideration  of  Stress  Singularity  and  Closing 
Pressure 


Fig.  A2.  Stress  Intensity  Factor  of  an  Infinite  Strip  Subjected  to 
a  Unit  Point  Load 


